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Abstract 

Abstract.  This  study  is  motivated  by  three-dimensional  ows  about  protrusions  and  cavities  with  an 
arbitrary  angle  between  the  external  ow  and  rigid  elements.  The  novel  type  of  a  \building  block”  cavity  ow 
isproposed  where  the  cavity  lid  moves  along  its  diagonal  (Case  A).  The  proposed  case  is  taken  as  a  typical 
representative  of  essentially  three-dimensional  highly  separated  vortical  ows  having  simple  single-  block 
rectangular  geometry  of  computational  domain.  Computational  results  are  compared  to  the  previous 
studies  where  the  lid  moves  parallel  to  the  cavity  side  walls  (Case  B).  These  3-D  lid-driven  cavity  ows  are 
studied  by  numerical  modeling  using  second-order  upwind  schemes  for  convective  terms.  The  volume  and 
plane  integrals  of  primary  and  transversal  momentum  are  introduced  to  compare  cases  in  a  quantitative 
way.  For  the  laminar  ow  in  the  cubic  cavity,  the  integral  momentum  of  the  secondary  ow  (which  is 
perpendicular  to  the  lid  direction)  is  about  an  order  of  magnitude  larger  than  that  in  Case  B.  In  Case  A,  the 
number  of  secondary  vortices  substantially  depends  on  the  Re  number.  The  secondary  vortices  in  the 
central  part  of  the  cavity  in  Case  A  distinguishes  it  from  Case  B,  where  only  corner  secondary  vortices 
appear.  For  a  rectangular  3-D3:  1  :  1  cavity  the  integral  momentum  of  the  secondary  ow  in  Case  A  is  an 
order  of  magnitude  larger  than  that  in  the  benchmark  cases.  The  ow  eld  in  Case  A  includes  a  curvilinear 
separation  line  and  non-  symmetrical  vortices  which  are  discussed  in  the  paper.  The  estimated  G.ortler 
number  is  approximately  4.5  times  larger  in  Case  A  than  that  in  Case  B  for  the  same  Re  number.  This 
indicates  that  in  Case  A  the  ow  becomes  unsteady  for  smaller  Re  numbers  than  in  Case  B.  For  developed 
turbulent  ow  in  the  cubic  cavity,  the  yaw  eect  on  ampli  cation  of  secondary  ow  is  as  strong  as  that  for  the 
laminar  ow  despite  the  more  complex  vortical  ow  pattern  in  benchmark  case  B. 
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THREE-DIMENSIONAL  FLOW  IN  CAVITY  AT  YAW 


ALEX  POVITSKY* 

Abstract.  This  study  is  motivated  by  three-dimensional  flows  about  protrusions  and  cavities  with  an 
arbitrary  angle  between  the  external  flow  and  rigid  elements.  The  novel  type  of  a  “building  block”  cavity 
flow  is  proposed  where  the  cavity  lid  moves  along  its  diagonal  (Case  A).  The  proposed  case  is  taken  as  a 
typical  representative  of  essentially  three-dimensional  highly  separated  vortical  flows  having  simple  single¬ 
block  rectangular  geometry  of  computational  domain.  Computational  results  are  compared  to  the  previous 
studies  where  the  lid  moves  parallel  to  the  cavity  side  walls  (Case  B).  These  3-D  lid-driven  cavity  flows  are 
studied  by  numerical  modeling  using  second-order  upwind  schemes  for  convective  terms.  The  volume  and 
plane  integrals  of  primary  and  transversal  momentum  are  introduced  to  compare  cases  in  a  quantitative  way. 
For  the  laminar  flow  in  the  cubic  cavity,  the  integral  momentum  of  the  secondary  flow  (which  is  perpendicular 
to  the  lid  direction)  is  about  an  order  of  magnitude  larger  than  that  in  Case  B.  In  Case  A,  the  number  of 
secondary  vortices  substantially  depends  on  the  Re  number.  The  secondary  vortices  in  the  central  part  of  the 
cavity  in  Case  A  distinguishes  it  from  Case  B,  where  only  corner  secondary  vortices  appear.  For  a  rectangular 
3-D  3:1:1  cavity  the  integral  momentum  of  the  secondary  flow  in  Case  A  is  an  order  of  magnitude  larger 
than  that  in  the  benchmark  cases.  The  flowfield  in  Case  A  includes  a  curvilinear  separation  line  and  non- 
symmetrical  vortices  which  are  discussed  in  the  paper.  The  estimated  Gortler  number  is  approximately  4.5 
times  larger  in  Case  A  than  that  in  Case  B  for  the  same  Re  number.  This  indicates  that  in  Case  A  the  flow 
becomes  unsteady  for  smaller  Re  numbers  than  in  Case  B.  For  developed  turbulent  flow  in  the  cubic  cavity, 
the  yaw  effect  on  amplification  of  secondary  flow  is  as  strong  as  that  for  the  laminar  flow  despite  the  more 
complex  vortical  flow  pattern  in  benchmark  case  B. 

Key  words,  three-dimensional  cavity  flow,  flow  at  yaw,  separated  flow,  incompressible  flow,  vortical 
flow,  laminar  and  turbulent  flow 

Subject  classification.  Fluid  Mechanics 

1.  Introduction.  The  lid-driven  cavity  flow  mimics  many  aeronautical,  environmental,  and  industrial 
flows  [10,  12,  2].  The  study  [2]  includes  a  detailed  description  of  cavity-type  flows  in  chemical  engineering. 
Authors  [10]  characterize  the  cavity  flow  problem  as  an  ideal  vehicle  for  studying  complex  environmental 
flows  in  a  simple  geometry.  Recent  aeronautical  motivation  comes  from  the  problem  of  airframe  noise. 
The  deployment  of  the  trailing-edge  flaps,  the  leading-edge  slats,  and  the  undercarriage  become  one  of  the 
important  sources  of  aircraft  noise  at  landing  phase  of  flight  [6,  20,  9].  Enclosures  of  quite  different  scales 
are  formed  by  wheel  wells,  brace  boxes,  and  pin  cavities  in  various  joints  linking  different  gear  components. 
Wheels  and  axles  comprise  a  series  of  short  cylinders  of  different  aspect  ratios  and  inclinations  relative  to  the 
flow  direction.  Therefore,  noise  suppression  associated  with  wheel-well  cavity  and  inter-wheel  gap  calls  for 
detailed  investigation  of  3-D  cavity-type  flows.  The  experimental  studies  [7,  8],  motivated  by  landing  gear 
noise,  demonstrate  that  acoustic  pressures  of  high  intensity  may  be  obtained  in  cavities  yawed  to  certain 
angles  even  if  the  flow  does  not  show  any  significant  oscillations  for  cavities  normal  or  parallel  to  the  mean 
flow  direction. 
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while  the  author  was  in  residence  at  ICASE,  NASA  Langley  Research  Center,  Hampton,  VA  23681-2199. 
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In  previous  studies,  the  cavity  flow  was  driven  by  a  lid  which  moves  in  the  direction  parallel  to  the  cavity 
side  walls.  This  flow  was  studied  numerically  [12,  10,  5],  experimentally  [13,  2],  and  has  become  a  benchmark 
case  for  developers  of  numerical  methods  [3,  22,  19],  to  name  a  few  in  each  category.  The  flow  is  almost  2-D 
at  low  Re  numbers,  transforms  to  the  3-D  flow  with  secondary  vortices  in  the  spanwise  plane  at  moderate  Re 
numbers,  and  bifurcates  from  a  stationary  flow  to  a  periodical  flow  at  a  critical  Re  number.  The  flow  loses 
its  symmetry  and  starts  to  induce  longitudinal  vortices  which,  in  turn,  lead  to  turbulence  [5,  10].  The  critical 
Re  number  depends  upon  depth-to- width  and  length-to- width  aspect  ratios.  This  problem  is  characterized 
by  quite  simple  geometry  and  boundary  conditions  and  complicated  vortical  structure  of  the  flow. 

In  real  life,  a  flow  is  likely  to  be  non-parallel  to  the  cavity  walls.  The  cavity  may  be  of  non-unity  aspect 
ratio  and  the  flow  character  may  range  from  steady  laminar  to  fully  developed  turbulent.  To  study  broader 
cases  than  that  described  above  and  keep  the  simple  geometrical  framework,  we  consider  the  lid  that  moves 
along  its  diagonal  (denoted  as  Case  A  here,  see  Fig.  la).  The  flow  impinges  in  the  spatial  angle  formed  by 
the  front  and  a  side  wall  that  leads  to  formation  of  the  system  of  vortices  caused  by  flow  separation.  In  spite 
of  the  lower  level  of  the  primary  flow  velocity  (in  the  direction  of  the  moving  lid),  the  numerous  secondary 
vortices  (in  the  perpendicular  direction)  appears  to  be  substantially  stronger  than  those  in  Case  B. 

We  investigate  this  flow  numerically  for  Re  numbers  corresponding  to  existence  of  steady  laminar  so¬ 
lution.  The  flowfield  and  integrals  of  primary  and  secondary  momentum  are  compared  to  those  in  the 
benchmark  Case  B.  The  location  of  vortices  and  their  dependence  on  the  Re  number  appears  to  be  quite  dis¬ 
tinctive  from  that  in  Case  B.  Next,  we  simulate  the  flow  in  the  1:1:3  depth-to- width-to-length  rectangular 
3-D  cavity  at  yaw.  The  yaw-related  enhancing  of  secondary  flow  leads  to  even  more  complex  pattern  of  sec¬ 
ondary  vortices  than  that  for  the  cubic  cavity.  Finally,  we  consider  the  developed  turbulent  cavity  flow  and 
obtain  the  amplification  of  the  secondary  flow  for  cavity  at  yaw.  Again,  the  yaw-related  amplification  of  the 
secondary  flow  leads  to  more  than  one  order  of  magnitude  increase  of  transversal  momentum  in  comparison 
to  Case  B. 

The  study  is  composed  as  follows.  In  Section  2,  the  problem  arrangement,  the  mathematical  model, 
boundary  conditions,  used  numerical  method,  and  validation  of  the  computer  code  are  described.  In  Section 
3,  the  features  of  primary  flow  are  discussed.  In  Section  4,  the  vortical  flow  patterns  of  secondary  flow  are 
presented  and  compared  to  those  in  Case  B.  In  Section  5,  the  study  has  been  expanded  to  fully  developed 
turbulent  flows  and  to  rectangular  (non-cubic)  cavities.  In  Section  6,  conclusions  and  future  work  are 
summarized. 

2.  Mathematical  model,  numerical  method,  and  code  validation.  The  governing  equations  are 
the  three-dimensional  steady  Navier-Stokes  equations.  The  boundary  conditions  describe  no-slipping  and 
non-penetrating  requirements  at  the  steady  walls  and  at  the  moving  lid  (see  Figs,  la,  b). 

In  Fig.  1,  the  considered  computational  volume  is  a  3-D  cube  with  its  edge  equal  to  unity.  In  Fig.  la 
(Case  A),  the  lid  moves  in  its  own  plane  in  the  direction  of  its  diagonal.  The  primary  flow  (in  the  direction 
of  moving  lid)  will  be  presented  in  the  main  diagonal  plane  CP  and  planes  CP2  and  CP 3,  which  are  parallel 
to  CP.  These  sections  are  2-D  rectangles  with  aspect  ratios  \[2  :  1,  1  :  1,  and  0.5  :  1,  respectively. 

To  calculate  the  Gortler  number,  the  velocity  profile  will  be  plotted  along  lines  ML  and  RL  in  the  CP 
plane  (see  Section  5).  The  former  line  crosses  the  boundary  layer  between  the  primary  vortex  and  the  cavity 
wall  while  the  latter  line  corresponds  to  the  radial  direction  of  the  primary  vortex.  The  flow  perpendicular 
to  the  direction  of  the  moving  lid  is  denoted  as  the  secondary  flow  and  will  be  shown  at  plane  PP  to  study 
flow  motion  from  the  cavity  bottom.  The  vortex  flow  pattern  will  be  shown  in  the  mid-plane  MP,  where 
y  =  0.5  (see  Section  4). 


2 


Fig.  1.  Cavity  flow:  a)  Case  A:  The  lid  moves  along  its  diagonal ;  b)  Case  B:  The  lid  moves  parallel  to  a  side  wall. 


For  Case  B  (Fig.  lb),  the  velocity  vector  field  in  planes  CP ,  PP,  and  MP  will  be  shown  for  comparison 
of  computational  results  to  Case  A.  These  planes  have  the  same  orientation  with  regard  to  the  direction 
of  the  moving  lid  as  those  in  Case  A.  For  code  validation  purposes,  the  velocity  profiles  along  the  vertical 
centerline  CL  for  various  numerical  grids  and  discretizations  of  governing  equations  will  be  discussed  later 
in  this  Section. 

The  governing  equations  are  discretized  on  uniform  and  non-uniform  three-dimensional  grids  using 
second-order  upwind  schemes  for  convective  terms  [11],  second-order  central  differences  for  diffusion  terms, 
and  second-order  approximation  for  pressure.  The  FLUENT  software  [1]  with  the  second-order  upwind 
discretizations  [4]  is  used  in  this  study.  FLUENT  uses  SIMPLE-type  [17]  algorithm  SIMPLEC  [21]  to  re¬ 
solve  velocity  and  pressure  coupling  and  algebraic  multigrid  (AMG)  to  solve  the  Poisson  pressure  correction 
equation. 

To  validate  the  code,  uniform  and  non-uniform  numerical  grids  are  exploited.  The  uniform  numerical 
grids  with  h  =  0.02  (51  x  51  x  51)  and  h  =  0.01  (101  x  101  x  101)  are  used.  The  non-uniform  grid  is  stretched 
near  rigid  boundaries.  The  smallest  grid  size  h  =  0.007,  the  ratio  of  grid  sizes  for  two  neighboring  nodes 
is  1.05.  The  non-uniform  part  of  the  grid  includes  20  rows  at  each  side;  therefore,  the  maximum  (the  20th) 
mesh  size  is  equal  to  0.177  and  the  non-uniform  part  occupies  2  x  23.14%  of  the  edge.  The  central  part  is 
covered  with  the  uniform  grid  (h  =  0.02).  The  total  number  of  grid  nodes  is  683,  which  is  3.2  times  less  than 
for  the  10 13  uniform  numerical  grid. 

For  Case  A,  the  profile  of  the  primary  velocity  (U  +  W)/y/2  along  the  line  CL  is  plotted  in  Fig.  2a.  The 
results  of  the  second-order  computations  are  closed  for  the  h  =  0.01  uniform  grid  and  for  the  non-uniform 
grid,  described  above.  Results  of  the  first-order  computations  on  the  same  grids  are  slightly  different  while 
0.3  <  y  <  0.6. 

For  Case  B,  the  computational  results  are  compared  to  the  pseudo-spectral  [14]  and  the  second-order 
[3]  numerical  solutions.  In  Fig.  2b,  the  computed  primary  velocity  W  along  the  centerline  is  presented  for 
Re  =  1000.  For  the  second-order  upwind  approximation,  the  obtained  velocity  profiles  are  closed  to  each 
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Fig.  2.  Validation  of  the  computer  code:  profile  of  primary  velocity  at  the  cavity  centerline  CL  :  a)  Case  A;  b)  Case  B. 

other  on  both  uniform  grids.  The  profiles  obtained  using  the  1013  uniform  grid  and  the  non-uniform  grid 
coincide  with  the  benchmark  results  [3,  14].  On  the  contrary,  the  first-order  computations  on  the  1013 
uniform  grid  show  discrepancy  with  the  benchmark  profile.  Computational  results  presented  in  the  next 
sections  were  obtained  in  the  non-uniform  grid  described  above. 

3.  Primary  flow.  Computed  velocity  vector  field  in  the  plane  CP  for  Case  A  is  presented  in  Fig.  3. 
This  section  represents  a  rectangle  with  the  aspect  ratio  y/2  :  1.  The  primary  vortex  is  shifted  towards  the 
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a. 


Fig.  3.  Case  A:  primary  flow  in  the  CPI  plane ,  Re  =  700 


Fig.  4.  Case  A:  primary  flow  in  cross-sections  with  various  aspect  ratio:  a)  CP 2  (aspect  ratio  1:1)  and  b)  CPS  (aspect 
ratio  0.5  :  1).  The  vector  scale  is  doubled  in  comparison  to  that  in  Fig.  3. 


upper  front  corner  of  the  cavity  in  comparison  to  the  2-D  flow  in  the  cavity  with  the  same  aspect  ratio. 

The  primary  flow  in  planes  CP2  and  CP 3  are  presented  in  Fig.  4.  In  plane  CP 2,  the  vortex  core  slowly 
rotates  at  the  cavity  central  part.  The  flow  current  is  directed  parallel  to  the  side  walls  and  along  the  cavity 
bottom.  In  plane  CP 3,  the  vortex  center  is  shifted  toward  the  upper  right  corner.  The  flow  is  directed 
upwards  in  the  most  of  the  section  and  does  not  contribute  to  the  primary  flow  (which  is  parallel  to  the  lid) . 
Recall  that  the  arrow  scale  is  doubled  in  Fig.  4  in  comparison  to  that  in  Fig.  3,  i.e.,  the  velocity  in  sections 
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Fig.  5.  The  integral  Is  of  primary  momentum 


CP2  and  CP 3  is  substantially  smaller  than  that  in  the  CP  plane. 

To  make  a  quantitative  comparison  of  the  primary  momentum  transfer  in  Cases  A  and  B,  the  volume 
and  surface  integrals  of  the  flow  momentum  in  the  direction  parallel  to  a  moving  lid  are  introduced.  Since 
the  lid  direction  in  Case  A  is  (i  +  k)  / y/2,  the  expression  for  local  momentum  of  primary  velocity  is  given 
by  =  0.5 {U  +  W )2.  In  Case  B,  the  momentum  of  primary  velocity  is  equal  to  M\\  =  W2.  The  volume 
integral  of  momentum  is  equal  to 


The  surface  integral  is  given  by: 


In  Fig.  5,  the  integral  Is  is  plotted  as  a  function  of  a  normalized  distance  n  measured  from  the  section 
CP.  In  order  to  plot  integrals  (2)  together  for  Cases  A  and  B,  the  distance  from  the  central  plane  is 
normalized  by  half-edge  1/2  in  Case  B  and  half-diagonal  y/2/2  in  Case  A.  In  both  cases,  the  Is  practically 
coincide  at  the  central  plane  (n  =  0).  In  Case  B,  the  Is  remains  constant  for  n  <  0.8  and  decreases  toward 
the  side  wall  for  larger  values  of  n.  In  other  words,  the  “boundary  layer”  thickness  is  equal  to  0.1  in  Case  B. 
In  Case  A,  the  Is  decreases  gradually  with  n.  The  Is  and  Iy  are  related  by: 

pn=  1 

(3)  Iv=  </>(n)Is{n)dn, 

Jn= 0 

where  the  weight  function  (j>{n)  =  2(1  —  n)  in  Case  A  and  (j>{n)  =  1  in  Case  B.  In  other  words,  the  variable 
section  area  in  Case  A  leads  to  a  non-trivial  weight  function  cj)(n)  and,  therefore,  the  Is  of  the  mean 
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Fig.  6.  Volume  integral  of  momentum:  a)Iy  and  b)  I ±. 

diagonal  section  CP  in  Case  A  has  doubled  weight  in  above  formulae.  For  both  cases,  the  volume  integrals 
of  momentum  (1)  are  presented  in  Table  1  and  in  Fig.  6.  Depending  on  the  Re  number,  the  integral  of 
momentum  is  20  —  30%  larger  in  Case  B  than  that  in  Case  A. 

The  length  of  a  2-D  CP  section  in  Case  A  is 

(4)  l  =  V2(l-n). 

For  instance,  n  =  0  for  CP;  n  =  1  —  1  / \/2  0.292  for  CP2;  and  n  =  1  —  0.5/\/2  0.646  for  CP3.  To 

compare  the  momentum  transfer  to  that  for  2-D  cavities  with  the  same  aspect  ratio,  the  computed  Is  for 
2-D  cavities  as  a  function  of  n  is  presented  in  Fig.  5.  The  2-D  cavity  hight  is  equal  to  unity  while  its  length 
is  computed  by  Eq.  (4).  In  the  2-D  case,  the  numerical  method  and  the  numerical  grid  used  are  the  same 
as  for  3-D  cases  (see  the  previous  Section).  The  integral  of  momentum  Is  for  the  2-D  \/2  :  1  cavity  is  135% 
higher  than  that  for  the  central  section  CP  of  Case  A.  The  difference  is  112%  for  CP2  and  29.9%  for  CP3. 
The  difference  is  caused  by  the  finite  span  in  3-D  cavities  and  variable  aspect  ratio  of  longitudinal  sections 
in  Case  A.  For  Case  B,  the  finiteness  of  the  cavity  span  leads  to  the  70%  difference  between  Is  for  CP 
section  and  the  2-D  square  cavity.  Note  that  for  the  rectangular  1:1:3  cavity  of  larger  span  (Section  5), 
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Fig.  7.  Flow  in  the  plane  perpendicular  to  the  direction  of  the  lid  (Re  =  700):  a)  Case  A;  the  vector  scale  is  doubled  in 
comparison  to  that  in  Fig.  3.  b)  Case  B;  the  vector  scale  is  four  times  larger  than  that  in  Fig.  3. 


the  difference  in  the  cross-section  reduces  to  21.5%. 

4.  Secondary  flow.  The  flow  in  the  plane  PP  is  shown  in  Figs.  7a, b  for  Cases  A  and  B,  respectively. 
In  Case  B,  the  vortices  appear  in  the  upper  and  lower  corners  of  the  section.  In  Case  A,  contrary  to  Case 
B,  a  vortex  pair  appears  near  the  symmetry  line.  At  the  symmetry  line  there  is  a  stagnation  point  located 
at  y  —  0.74.  From  the  stagnation  point,  the  flow  is  directed  up  towards  the  lid  and  down  towards  the  cavity 
bottom.  The  vortex  size  and  strength  are  considerably  larger  than  those  in  Case  B.  (The  vector  scale  is 
doubled  in  Fig.  7b  in  comparison  to  that  in  Fig.  7a)). 

To  make  a  quantitative  comparison  of  transverse  momentum  transfer  in  Cases  A  and  B,  the  momentum 
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Table  1 

Volume  integrals  of  momentum  for  3-D  cubical  cavity 


Re  number 

Case  A 

Case  B 

Iv 

i± 

Iv 

i± 

100 

0.34E-01 

0.17E-02 

0.358E-01 

0.14E-03 

400 

0.25E-01 

0.27E-02 

0.307E-01 

0.436E-03 

700 

0.216E-01 

0.283E-02 

0.270E-01 

0.389E-03 

1000 

0.189E-01 

0.253E-02 

0.247E-01 

0.329E-03 

2000 

0.163E-01 

0.239E-02 

0.183E-01 

0.175E-03 

in  the  direction  perpendicular  to  a  moving  lid  were  considered.  Since  the  transverse  direction  in  Case  A 
is  (i  —  k)/y/2,  the  expression  for  local  transverse  momentum  is  given  by  M±  =  0.5 (U  —  W)2 .  In  Case  B, 
the  transverse  momentum  is  equal  to  M±  =  U2.  The  volume-averaged  value  of  integral  of  momentum  is 
computed  by 

(5)  W/v"- 

Integral  I±  is  about  one  order  of  magnitude  larger  in  Case  A  than  that  in  Case  B  (Fig.  6b).  For 
Re  =  2000,  the  I±  is  of  14.66%  of  Iy  in  Case  A  whereas  the  I±  is  of  only  0.95%  of  Iy  in  Case  B.  In  fact, 
much  stronger  transverse  flow  motion  in  Case  A  confirms  this  integral  evaluation  (compare  Figs.  7a  and  7b). 

To  study  the  transverse  flow  motion  in  the  middle  plane  MP,  the  velocity  vectors  are  shown  in  Figs.  8a-d 
and  Figs.  9a-d  for  Cases  A  and  B,  respectively,  and  the  range  of  Reynolds  numbers  within  100  <  Re  <  2000. 
The  plotted  vector  scale  is  the  same  in  all  Figs.  8-9  and  equal  to  this  in  Fig.  7b.  For  Re  =  100,  the  flow  is 
parallel  to  the  lid  direction  and  no  visible  vortices  are  observed  in  Cases  A  and  B.  For  Re  =  400,  the  vortex 
couple  appears  in  the  MP  corners  in  both  Cases  A  and  B  and  these  vortices  are  designated  as  1  in  Fig. 
8b.  In  Case  B,  the  only  vortex  couple  is  visible  for  larger  Re  numbers  (Figs.  9c, d).  These  vortices  elongate 
towards  the  symmetry  line  with  increase  of  the  Reynolds  number.  As  opposed  to  Case  B,  in  Case  A  the 
number  of  vortex  pairs  increases  with  the  Re  number  and  five  couples  of  vortices  are  visible  and  designated 
as  1,2,  3, 4,  and  5  in  Fig.  8d  {Re  =  2000). 

5.  Rectangular  cavity  and  turbulent  flow.  The  computational  results  for  a  3-D  rectangular  cavity 
with  1:1:3  depth-to-width-to-length  ratio  is  considered.  The  cavity  lid  is  1:3  rectangle;  the  Re  number  is 
equal  to  700.  Three  cases  are  modeled:  (i)  the  lid  moves  along  its  diagonal  (Case  A);  (ii)  the  lid  moves  along 
its  long  side  (Case  Bl);  and  (iii)  the  lid  moves  along  its  short  size  (Case  B2). 

In  Case  Bl,  the  strong  vortices  appear  at  the  front  corners  where  the  flow  impinges  into  the  front  wall 
of  the  cavity.  The  rest  of  computational  volume  includes  weak  vortices  near  side  walls  and  strong  back  flow 
along  the  cavity  axis.  In  Case  B2,  the  vortices  occupy  the  space  near  side  walls  up  to  the  back  wall.  In  Case 
Bl,  the  volume-averaged  integral  Iy  is  larger  and  I±  is  smaller  than  those  in  Case  B2  (see  Table  2).  In  Case 
A,  the  curvilinear  stagnation  line  initially  coincides  with  the  diagonal  of  PM  and  then  deviates  towards 
the  side  wall.  Non-symmetrical  vortices  with  respect  to  the  diagonal  of  the  lid  appear  at  both  sides  of  the 
diagonal  in  the  back  part  of  the  cavity.  The  value  of  integral  of  primary  momentum  Iy  in  Case  A  is  between 
those  for  cases  Bl  and  B2.  On  the  contrary,  the  value  of  integral  of  secondary  momentum  I±  in  Case  A  is 
approximately  one  order  of  magnitude  larger  than  those  in  Cases  Bl  and  B2. 


9 


Table  2 

Volume  integrals  of  momentum  for  rectangular  cavity,  Re=700 


Cas 

Iv 

e  A 

I± 

Cas< 

Iv 

B1 

I± 

Cas< 

Iv 

B2 

I± 

0.403E-01 

0.324E-02 

0.460E-01 

0.217E-03 

0.359E-01 

0.411E-03 

The  goal  of  this  study  is  to  investigate  the  cavity  flow  for  Re  numbers  corresponding  to  existence  of 
steady  solution  for  such  a  flow.  Nevertheless,  the  Gortler  numbers  for  Cases  A  and  B  based  on  steady 
flowfield  are  computed  to  indicate  that  the  transition  to  unsteady  flow  may  occur  for  smaller  Re  numbers 
in  Case  A. 

The  profile  of  vertical  velocity  V  along  the  middle  line  ML  is  shown  in  Fig.  12  for  Cases  A  and  B 
for  Re  =  700.  In  Case  A,  the  abscissa  has  been  scaled  by  y/2.  The  absolute  value  of  maximum  velocity  is 
0.65  in  Case  A  and  0.43  in  the  Case  B.  Along  the  line  RL ,  the  ratio  of  the  maximum  absolute  velocities 
is  approximately  equal  to  1.5  as  well.  The  thickness  of  boundary  layer  estimated  by  the  distance  from  the 
wall  to  the  point  of  maximum  velocity  is  the  same  in  Cases  A  and  B.  However,  the  actual  thickness  of  the 
boundary  layer  in  Case  A  is  y/2  larger  because  of  scaling.  The  Gortler  number  is  based  on  the  boundary 
layer  Reynolds  number,  the  boundary  layer  thickness  and  the  radius  of  curvature  of  the  wall  [10]: 


(6) 


G  = 


Re25 

rc 


where  Re  =  UmaxS/v  is  the  boundary  layer  Reynolds  number,  S  is  the  boundary  layer  thickness,  and  rc  is 
the  radius  of  curvature  of  the  wall.  Authors [10]  assume  that  rc  is  equal  to  the  half-length  of  cross-section, 
i.e. ,  rc  =  1/2  in  Case  B  and  rc  =  y/2/2  in  Case  A.  The  ratio  of  Gortler  numbers  in  Cases  A  and  B  is  given  by: 
Ga/Gb  =  1.52-\/2  /y/2  =  4.5.  Therefore,  the  unsteady  Gortler  vortices  may  appear  at  smaller  Re  numbers 
in  comparison  with  Case  B  and  may  convect  to  the  peripheral  parts  of  the  cavity  by  transverse  flow. 

For  the  turbulent  flow  in  cubical  cavity  the  Navier-Stokes  equations  are  replaced  with  the  Reynolds- 
averaged  steady  equations  and  the  Reynolds  stress  model  of  turbulence  [15].  The  Launder  and  Spalding  [16] 
wall  functions  are  used  for  boundary  conditions  for  the  Reynolds  stresses. 

For  Re  >  2000,  the  cavity  flow  becomes  unsteady.  Previous  experimental  [18]  and  computational  [23] 
studies  of  Case  B  for  3000  <  Re  <  10000  show  that  the  cavity  flow  is  unsteady  with  high-  (turbulent) 
and  low-frequency  (organized  structure)  fluctuations  of  velocity  around  time- averaged  values.  While  for 
Re  «  3000  the  fluctuations  are  due  to  low-frequency  meandering  of  Gortler  vortices,  the  break-down  of 
organized  Gortler  vortices  is  observed  for  Re  =  10000  and  the  high-frequency  turbulent  fluctuations  become 
responsible  for  momentum  transfer.  This  type  of  unsteadiness  is  modeled  using  the  RANS  approach  in  the 
current  study.  Still,  the  presence  of  single  vortices  sweeping  through  the  computational  volume  is  quite 
possible  and  may  affect  the  momentum  transfer  [18].  An  assumption  about  the  presence  of  steady  solution 
to  RANS  equations  is  adopted  in  this  study  and  the  effect  of  low-frequency  unsteady  vortex  shedding  will 
be  considered  in  our  future  research. 

The  flowfield  in  the  cross-section  CP  is  shown  in  Fig.  13  for  Re  =  40000.  For  Re  =  10000,  the  flowfield 
shows  similar  vortical  pattern.  To  compare  with  the  laminar  flowfield  in  Case  A,  each  of  two  central  vortices 
is  split  into  three  unequal  vortices  (compare  Figs.  7a  and  13a).  The  corner  flow  separation  leads  to  formation 
of  vortices  in  four  corners  of  section  CP ,  which  are  not  visible  in  the  laminar  case.  In  Case  B,  the  corner 
vortices  are  split  and  the  vortical  area  is  shifted  toward  the  vertical  centerline  (compare  Figs.  7b  and  13b). 
For  Case  A,  five  pairs  of  vortices  are  observed.  Their  location  in  the  mid-section  MP  is  similar  to  that  for 
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Table  3 

Volume  integrals  of  momentum  for  turbulent  cubical  cavity  flow 


Re 

Cas 

Iv 

e  A 

I± 

Cas 

Iv 

le  B 

I± 

10000 

40000 

0.717E-02 

0.365E-02 

0.120E-02 

0.681E-03 

0.967E-02 

0.414E-02 

0.162E-03 

0.540E-04 

laminar  flow  (Re  =  2000)  and  turbulent  flow  (compare  Figs.  8d  with  14a).  In  Case  A  and  Re  =  40000,  the 
integral  of  momentum  of  transverse  motion  makes  18.6%  of  the  integral  of  primary  momentum  (see  Table 
3).  For  Case  B,  the  flowfield  pattern  is  quite  different  for  turbulent  and  laminar  flows  (compare  Figs.  9  and 
14b).  Numerous  vortices  appear  at  the  front  and  back  parts  of  the  section  MP.  Nevertheless,  there  is  no 
significant  increase  of  transversal  part  of  momentum  as  a  fraction  of  primary  momentum  since  the  integral 
of  momentum  of  transverse  motion  I±  is  of  1.3%  of  the  integral  of  primary  momentum  Iy.  For  Re  =  10000, 
the  ratio  between  different  components  of  flow  momentum  remains  approximately  the  same  (see  Table  3) 
while  the  values  of  Iy  and  I±  are  doubled  in  comparison  to  those  for  Re  =  40000. 

6.  Conclusion  and  future  plans.  The  cavity  flow  driven  by  the  cavity  lid  moving  along  its  diago¬ 
nal  is  simulated  using  the  the  second-order  upwind  discretizations  of  convection  terms.  This  flow  mimics 
environmental  and  aeronautical  flows  where  the  governing  flow  direction  is  not  parallel  to  cavity  side  walls. 
The  considered  flow  (Case  A)  is  substantially  different  from  that  in  the  benchmark  case  where  the  cavity 
lid  moves  parallel  to  the  side  walls  of  the  cavity  (Case  B).  Enhanced  transverse  fluid  motion  in  the  direc¬ 
tion  perpendicular  to  the  direction  of  moving  lid  is  observed  in  Case  A.  The  integral  of  momentum  in  this 
direction  is  about  one  order  of  magnitude  higher  than  that  in  Case  B.  The  number  and  location  of  vortices 
in  the  mid-plane  MP  the  lid)  is  quite  different  from  that  in  Case  B,  where  only  corner  vortices  appear.  In 
Case  A,  the  number  of  vortices  increases  with  the  Re  number  and  these  vortices  appear  at  the  central  part 
of  the  mid-plane. 

The  flow  in  the  3-D  elongated  cavity  driven  by  the  cavity  lid  has  the  curvilinear  separation  line  and 
non-symmetric  vortices  in  the  mid-plane.  Again,  the  integral  of  transverse  momentum  is  approximately  one 
order  of  magnitude  higher  than  those  in  benchmark  cases  B,  where  the  lid  moves  parallel  to  either  long  or 
short  side  walls. 

For  diagonal-driven  fully  developed  turbulent  cavity  flow  in  a  cube,  modeled  by  the  Reynolds  Stress 
model  of  turbulence,  the  vortical  pattern  is  qualitatively  similar  to  that  for  the  laminar  flow  at  Re  =  2000. 
For  Case  B,  significantly  more  secondary  vortices  are  observed  in  planes  PP  and  MPa  than  for  laminar  flow 
case.  Nevertheless,  the  transverse  momentum  for  Case  A  remains  more  than  one  order  of  magnitude  higher 
than  that  for  Case  B,  i.e. ,  the  strong  yaw-related  amplification  of  transversed  flow  holds  for  turbulent  Re 
numbers. 

For  the  main  diagonal  plane  in  Case  A,  the  Gortler  number  appears  to  be  4.5  times  larger  than  that 
for  Case  B  for  the  same  Re  number.  Therefore,  the  unsteady  Gortler  vortices  may  appear  at  smaller  Re 
numbers  than  in  Case  B.  Although  the  second-order  discretizations  is  accurate  enough  for  the  viscous  and 
steady  flow,  higher  order  numerical  schemes  are  required  for  wavy  and  unsteady  transition  phenomena.  The 
high-order  numerical  modeling  of  transition  to  unsteady  flow  including  Gortler  vortices  will  be  conducted  in 
our  future  research. 
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Fig.  8.  Flow  in  3-D  cubical  cavity  in  Case  A:  a)  Re  =  100,  b)  Re  =  400,  c)  Re  =  700,  and  d)  Re  =  2000. 
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Fig.  9.  Flow  in  3-D  cubical  cavity  in  Case  B:  a)  Re  =  100,  b)  Re  =  400,  c)  Re  =  700,  and  d)  Re  =  2000. 


Fig.  12.  Velocity  profile  along  the  line  ML. 
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Fig.  13.  Turbulent  cavity  flow  in  the  plane  perpendicular  to  the  direction  of  the  lid.  a)  Case  A:  the  vector  scale  is  six 
times  larger  than  that  in  Fig.  3;  b)  Case  B:  The  vector  scale  is  twelve  times  larger  than  that  in  Fig.  3. 


